Influence of vortices on the magnetic resonance in cuprate superconductors 
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We investigate several theoretical possibilities for the suppression in a c-axis magnetic field of the magnetic 
resonance recently observed in inelastic neutron scattering experiments on YBa2Cu306.e- We find that 
neither the Doppler shift of the quasiparticle states caused by supercurrents outside the vortex core, nor 
an assumed spatially uniform suppression of the coherence factors or spectral gap due to the applied field, 
can account for the observed effect. In contrast, suppressing the gap or the coherence factors in the vortex 
core to zero is consistent with the data. We demonstrate that an even simpler description of the data can 
be achieved by assuming that the resonance is not supported within an effective radius £ e // around each 
vortex, where £ e f / is the sum of the superconducting and spin-spin correlation lengths. We use this simple 
idea to predict the doping dependence of the field suppression. 
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One of the more intriguing developments in the field 
of high temperature cuprate superconductivity has been 
the observation by inelastic neutron scattering (INS) ex- 
periments of a siiarp magnetic resonance in the super- 
conducting statetl. Recently, it was found that a c-a: 
magnetic field suppressed the intensity of this resanano 
as predicted from an analysis of specific heat datatlJSmce 
the same effect was not observed for in-plane fieldsQ, this 
indicates that the resonance is sensitive to the presence 
of Abrikosov vortices, and thus intimately connected to 
the nature of the superconducting ground state. This 
has obvious implications for microscopic theories of the 
resonance. 

In this paper, we consider a model where the resonance 
is treated as a particle-hole bound state in a d-wave su- 
perconductor, with calculations performed within linear 
response theory (RPA) . Several effects of the vortices are 
considered. First, we calculate the influence of the super- 
currents circulating around the vortices on the resonance. 
We find that this only leads to a broadening of the reso- 
nance in energy; the integrated weight remains the same, 
in conflict with experiment. Second, we study the effect 
of a spatially uniform suppression of the (Ak Ak+cj) cor- 
relator which enters the coherence factors of the spin sus- 
ceptibility (where Q is the antiferromagnetic wave vec- 
tor at which the resonance is peaked). Such a suppres- 
sion is speculative, but could be a result of dephasing of 
the pairing in a c-axis field due to the vortices, as ob- 
served in Josephson plasma resonance experimentsu. We 
find that although this does lead to a suppression of the 
integrated weight as observed experimentally, the effect 
causes the resonance to shift to higher energy, in conflict 
with experiment. Third, an assumed (field induced) spa- 
tially uniform suppression of the gap magnitude causes 
the resonance to shift to lower energy, also in conflict 
with experiment. 

This leads us to consider the effect of the vortex cores 
themselves. We observe that if the resonance is not sup- 
ported in the vortex cores, then the resulting field depen- 



dence is in reasonable agreement with experiment. We 
consider three possibilities for the suppression of the res- 
onance in the vortex core regions: (a) the suppression of 
the gap magnitude in the core, (b) the suppression of the 
(AA) correlator in the core, and (c) the absence of quasi- 
particles in the core. Any of these three possibilities give 
a good account of the data. We use case (c) to estimate 
the doping dependence of the field suppression effect. 

To calculate the influence of the supercurrents around 
the vortices on the resonance in the spin-spin correlation 
function, we approximate the superflow by a circular flow 
around the vortex center. The corresponding local su- 
permomentum, p s , is proportional to the gradient of the 
phase, p s = fte^/2r. This is a good approximation for 
the experiments considered here, where the intervortex 
spacing is smaller than the penetration depth and large 
compared to the coherence length. 

In the intervortex regions, the variation of the order pa- 
rameter and of the superflow occurs on a scale large com- 
pared to the spin-spin correlation length, which amounts 
to only a few lattice constants as determined from the 
momentum width of the resonance. Consequently, we 
determine the RPA susceptibility in the intervortex re- 
gion at each point of the unit cell of the vortex lattice in 
the presence of the local superflow, 
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The bare susceptibility Xo(k>) Q> Ps)i is determined as 

xo(^q, Ps ) = -E e +K _ K +tT 
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where the excitation spectrum in the presence of a su- 
perflow with momentum p s is given by 
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with = (f k+p , - £k-p s )/2, £ k - (£k+ Ps + £k-p s )/2, 
and 

^k = I ± ^k = ± + Ak (4) 

as coherence factors. The factor a (1 in the current 
case) will be discussed later. The spatial average of 
the Doppler shifted susceptibility over the intervortex 
region of the vortex lattice unit cell is then calculated, 
x(co,Q) =< x(w,Q,p s (R)) >r. We evaluated Eq. (§) 
for a 512 x 512 grid of k-points and pcrfomed the spatial 
average over 320 R-points. Note we use a normalization 
equal to the intervortex area, not the total area. The 
contribution from the vortex cores will be discussed later 
in the paper. 

Calculations were performed using a model quasiparti- 
cle dispersion in the superconducting state motivated by 
photoemission measurements!! Similar dispersions were 
found to give a good description of the zero field INS 
data, including the incomniensurate structure observed 
at energies below resonance!! A d-wave superconducting 
gap proportional to cos(k x a) ~ cos(k y a) was assumed, 
with a maximum value of A =i-29 meV as determined 
from recent STM measurements!! A broadening factor, 
r, of 2 meV was employed, and a temperature of 13 K. 

In Fig. [l] we show our results for the effect of the cir- 
culating supercurrents on the resonance. The exchange 
coupling, Jq, is fixed to give a resonance at 34 meV for 
zero magnetic field. For the spatial average we assumed 
a lower cut-off at the value ^ — 2a (vortex core radius) 
and an upper cut-off at the value R = 25a (radius for en- 
closing one flux quantum at 7 T), where a is the Cu-Cu 
distance. The results are insensitive to the lower cutoff. 
As Fig. |l| shows, the supercurrent has three effects: (a) 
it shifts the position of the resonance to slightly lower 
energy, (b) it broadens the resonance, and (c) it reduces 
the magnitude of the resonance at the peak energy. Also 
shown are the energy integrated susceptibilities, which 
demonstrate that the integrated weight between and 
pa 2A is conserved. These findings are in apparent con- 
tradiction with the experimental facts, which are that 
the resonance does not shift, nor broaden, and that the 
integrated weight is reduced by about 15% at 7 T [Ref. 

We have also tested a number of other dispersionsB, 
and a variety of assumed values for A and J. Although 
the amount of broadening is somewhat sensitive to these 
details, we find that the integrated weight is always ap- 
proximately conserved. An example is given in the right 
panel of Fig. |l|, where we find virtually no effect of the 
Doppler shift on the susceptibility. 

We also checked if an assumed field induced (spatially 
uniform) reduction of the gap magnitude accounts for the 
observed effect. Our result is shown as the dotted line in 
Fig. H compared to the zero field result (full line). The 
integrated weight is suppressed in this case (the left panel 
of Fig. shows the reduction of the integrated weight 
versus A^(iT)/A 2 (0)). To obtain the observed 15% 
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FIG. 1. Influence of Doppler shifts due to supercurrents 
on Im x(w,Q). The full lines correspond to zero magnetic 
field, the dashed ones to a field of 7 T. The thin lines are 
the energy integrated susceptibilities (scaled by a constant so 
that they could be plotted on the same graph). Although the 
magnitude of the resonance at the peak energy is suppressed, 
its integrated weight is not. The parameters employed were 
A=29 meV, T=2 meV, T=13 K. In the left panel we use J = 
357 meV, and the dispersion taken from Ref. In the right 
panel we use J = 142 meV and the dispersion tb2 taken from 
Ref. ^| for this case, the Doppler shift has virtually no effect. 

reduction in weight at 7 T would require reducing the gap 
from 29 to 20 meV. This reduction is substantially larger 
than would be indicated by the upper critical field (45 T), 
and the reduction appears to have the wrong functional 
dependence on H. Moreover, this gap reduction shifts the 
resonance to considerably lower energy, in contradiction 
with experiment. 

As a third mechanism, we studied a (spatially uniform) 
suppression of the (AkAk+cj) correlator in the CkCk+Q 
coherence factors (by reducing a to less than 1 in Eq. @). 
The motivation for this is that phase fluctuations induced 
by the vortices are known to lead to a dephasing of the 
layers, and the resonance will be sensitive to this since 
it involves c-axis coupling (it is peaked at k z = Tt/d, 
where d is the separation of nearest neighbor CuO lay- 
ers). The observed decoupling inferred from the field 
dependence of the Josephson plasmonQ, though, is prob- 
ably due to the weaker bilayer-bilayer coupling, which 
is also consistent with small mesa experiments!!!! There- 
fore, at the current time, it is not known whether the two 
layers within a bilayer are dephased or not (though this 
could be determined from the field dependence of c-axis 
infrared conductivity measurements, where a feature is 
seen attributed to an optical Josephson plasmonO). For 
now, though, we will .assume that this is a possibility, and 
test its consequencesE! 

In Fig. |H we compare the zero field result to the 
same result, but with the correlator reduced by 15% 
(a = 0.85). This leads to a large reduction of the in- 
tegrated weight, as seen experimentally (in Fig. ||, we 
plot the integrated weight versus a). We note that the 
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FIG. 2. Comparison of zero field susceptibility with the 
same susceptibility, but with reduced gap magnitude, A. Dot- 
ted line: assuming a spatially uniform reduction from 29 meV 
to 20 meV; dashed line: assuming a reduction to zero in 24% 
of the vortex unit cell area representing the cores (the uniform 
zero gap response is shown in the inset). In both cases the 
weight of the resonance is reduced, but for the spatially uni- 
form case, the resonance is shifted considerably downwards in 
energy. The parameters used are the same as in Fig. |lj 

experimental suppression goes like 1 — H/H*, where H* 
is a Jiumber not much lower than H C 2, the upper critical 
fieldo. Based on quantum Ginzburg-Landau theory, a re- 
duction of the (A A) correlator proportional to \ — H/H C 2 
is expected. Therefore, it is reasonable to suppose that 
the relative experimental suppression goes like a. This 
suppression is in good agreement with the calculation, 
as can be seen in Fig. ||. We note, however, that the 
position of the resonance shifts to higher energies, in dis- 
agreement with the data. It would be coincidental if this 
energy shift was exactly canceled by an assumed shift of 
the superconducting gap to lower energies by the field. 

Let us now consider the effect of the vortex cores. The 
fact that the experimental suppression goes like 1—H/H* 
is highly suggestive of a vortex core effect, as originally 
noted by Dai et alu. Therefore, we assume that the 
resonance is not supported in the region of the vortex 
core. This assumption is based on five facts: a) the con- 
siderable momentum width of the resonance shows that 
the corresponding spin excitations have a decay length 
of only two lattice constants, which is smaller than the 
coherence length; thus the resonance will be sensitive to 
variations of the order parameter on the coherence length 
scale; b) the resonance at zero field only exists in the su- 
perconducting state, and disappears in the normal state; 
c) coherence peaks in the single particle density of states 
at the gap edge wece not found in the core region in 
STM measurements^; this would modify the 2A-edge in 
Xo (Eq. ||) and suppress the resonance; d) in underdoped 
materials, missing subgap states point towards a loss of 
quasiparticle weight due to a pseudogap in the vortex 
core;E^I e) the dip feature in the tunneling density 
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FIG. 3. Comparison of zero field susceptibility with the 
same susceptibility, but with the (AA) correlator in the nu- 
merator of Eq. 2 reduced (a < 1). Dotted line: assuming a 
spatially uniform reduction by 15%; dashed line: assuming a 
reduction to zero in 20% of the vortex unit cell area repre- 
senting the cores (the uniform zero a response is shown in the 
inset). In both cases the weight of the resonance is reduced, 
but for the spatially uniform case, the resonance is shifted 
considerably upwards in energy. The parameters used are the 
same as in Fig. |l| 

of states, thought to be due to the coupling of quasiparti- 
cles tq-the resonanced, is not observed in the vortex core 
regiontil 

In Figs. H and |[ we show in the insets the suscepti- 
bility for zero A, and for zero (AA) correlator. In both 
cases, the resonance is strongly suppressed. In the main 
panels, we show as dashed curves the results for the case 
when we use the curves in the insets for the vortex core 
regions, and the full curves (zero field results) for the in- 
tervortex regions. The latter is justified since we found 
above that the Doppler shift has a negligible effect on the 
integrated intensity. In both cases, the resulting curves, 
calculated for a 15% reduction in total integrated weight, 
reproduce very well the experimental finding of no shift 
or broadening of the resonance. One problem is that 
the observed weight suppression would require a vortex 
core region filling 24% of the total area for the first case, 
and 20% for the second case. Both values are somewhat 
larger than what one would expect at 7 T, as discussed 
below, especially for the first case. 

There is a way of testing whether either of these two 
scenarios is correct. In the case where the gap is re- 
duced in the cores, then extra weight would show up at 
lower energies, whereas for the case where the correlator 
is reduced in the cores, extra weight would show up at 
higher energies. Although our calculations indicate that 
the corresponding changes between and 7 T are perhaps 
within experimental error bars, we encourage experimen- 
talists to look for extra weight, both in the region around 
20 meV and around 50 meV. 
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FIG. 4. Energy integrated weight (from uj = to ui = 50 
meV, normalized to the zero field value) from Figs. ^ and || 
Left: versus A 2 (i?)/A 2 (0), where A is the maximal d-wave 
gap magnitude. Right: versus a, where a is the prefactor 
in front of the (AA) correlator in the numerator of Eq. 2. 
The dotted lines are the expected behavior if the normalized 
weight was equal to (1 — H/ H C 2). 



to the core region (this was a motivation for the a < 1 
calculations). We suggest that this may lead to a de- 
struction of quasiparticle excitations in the vortex core 
region similar to what happens in the pseudogap state. 
The absence of quasiparticle peaks as well as the neu- 
tron resonance in the core region is consistent with the 
notion that both these spectral features require substan- 
tial local phase correlations^. While this conjecture is at 
this stage admittedly speculative, we believe it deserves 
further experimental and theoretical investigation. 

This work was supported by the U.S. Dept. of Energy, 
Office of Science, under Contract No. W-31-109-ENG-38. 



We note that our formalism assumes the presence of 
quasiparticle states, and thus we do not expect a reso- 
nance to exist if quasiparticles do not exist. A simpler 
idea to those discussed in the previous two paragraphs is 
to simply assume that the weight is zero in the considered 
energy range, consistent with a spin gap of ~ 50 meV in- 
side the cores. From this, we can estimate the effective 
radius of the core, £ e / / , from the ratio of the INS weight 
to that at zero field. Since the ratio of the vortex core 
area to the total area is equal to Bir^ff/^o, to obtain 
an effect of 15% at 7 T in underdoped YE^CusOy-^, 
we need a £ e / / of 37-38 A. This is somewhat larger than 
the estimated coherence length of 27 A (H c2 of 45 T), 
but it should be noted that £ e // should approximately 
be the sum of the superconducting and magnetic cor- 
relation lengths. The latter is approximately 10 A (as 
measured by the momentum width of the resonance), so 
the two estimates are in agreement. 

For optimally doped compounds, both coherence 
lengths are shorter. Assuming £ e ff «20-24 A, our pre- 
diction for this case is that the suppression of the total 
weight at 7 T will only be 4-6%. Going further to the 
overdoped regime, the superconducting coherence length 
increases, leading to an increase of the sensitivity of the 
resonance with magnetic field again. Further underdop- 
ing, though, should lead to an even more dramatic reduc- 
tion, as both coherence lengths are expected to increase 
as the doping is reduced. In fact, we would argue that 
the field dependence of the resonance at various dopings 
would be a good measure of the doping dependence of 
the superconducting coherence length, and it would be of 
great interest to correlate STM and INS measurements 
on the same samples. 

We wish to conclude this paper with the fol- 
lowing speculation, motivated by the above re- 
sults. As documented by angle-resolved photoemission 
measurement*!-], quasiparticle-like peaks in the spectral 
functions are present only below T c , the onset temper- 
ature of phase coherence. The superconducting phase 
is singular at the vortex core, and therefore the phase 
correlations are strongly suppressed between points close 
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